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We non-perturbatively calculate the renormalization factor of the static axial vector current in 0(a) improved 
quenched lattice QCD. Its scale dependence is mapped out in the Schrodinger functional scheme by means of a 
recursive finite-size scaling technique, taking the continuum limit in each step. We also obtain Z^''* for Wilson 
fermions in order to renormalize existing unimproved data on F^'^^ non-perturbatively. 



1. INTRODUCTION 

The decay constant Fb, governing the lep- 
tonic decay of a B-meson, is a quantity of much 
phenomcnological interest. However, as heavy 
flavours on the lattice escape a direct numerical 
treatment, the static approximation represents a 
valuable alternative, since this effective theory 
has simplified dynamics and describes the large 
mass limit of the theory. Yet the problems of this 
framework have been twofold in the past. Owing 
to the infinitely heavy b-quark, (i) the renormal- 
ization properties of the static theory are differ- 
ent, i.e. the renormalization constant Z^^^ of the 

axial current in (^r''*)o = ^a'^Hm) "0(17075 V'b''* 
becomes scale (^) dependent, thereby entailing 
an additional uncertainty, and (ii) MC computa- 
tions of the matrix element itself are difficult. 

Here, we solve (i) by matching through the 
renormalization group invariant (RGI) operator: 
The non-perturbatively calculated //-dependence 
of any renormalized matrix element of Aq**** , $ = 
^stat ^ (f |(Af|;''*)o| i), together with the value of 
Z^'^* at some (low) /i, yields a factor leading to 
the RGI $. It is then related to the 'matching' 
scheme, which approximates the relativistic the- 
ory up to 1/m-corrections. For a more detailed 
explanation of the overall strategy we refer to [Q 
and for a full report on our work to H . 



*Based on a talk presented by J.H. at the conference LAT- 
TICE '02, June 24 - 29, 2002, MIT, Cambridge MA, USA. 
tSupported by the EU under HPRN-CT-2000-00145. 



2. STRATEGY 



Our determination of ZJ^*^* and its running 
with /i uses the Schrodinger functional (SF) j|] 
as intermediate scheme Q . The procedure resem- 
bles alpha's quark mass renormalization 

1. Choose a proper combination X of SF corre- 
lation functions such that, with ATr = Z'^^^X 



finite and X = + X^^^^g + • • • , the condi- 
tion Xji = X^^^f non-perturbatively defines Z^^'^*, 
running with the SF's box size L = l//i. 

2. Map out the L-dependence recursively by em- 
ploying the finite-size step scaling function (SSF) 
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u = gl^{L), (1) 



computed on the lattice and extrapolated to the 
continuum as illustrated in fig. |l|. 

3. Evolve Z^p^^'^ from initially large L (low fx) to 
small L (high ^) by repeatedly applying [c^^'^*]"^; 
from there, continue with the perturbative high- 
energy behaviour to arrive at (scale and scheme 
independent) RGI matrix elements 
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$RGI = *(Ai) X {2hQ-f) 
X exp I - y dg 

of (A||;''*)o, with g = 5sf(-^ — 1/m) s-^^d universal 
coefficients 60 = ll/(47r)^ and 70 = -l/(47r^). 
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Finally, for sufficiently high the RG invariant 
^RGi can be converted to the 'matching' scheme 
via a perturbative factor $match(M)/*J'RGi- 
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Figure 1. Lattice step scaling function and its 
continuum extrapolation for some selected u. 

3. SURVEY OF THE RESULTS 

The simulation and analysis as well as the SF 
setup are basically analogous to except that 
the boundary coefficient ct is set to 2-loop ||]. 
While the inclusion of static quarks and their im- 
pact on 0(o) improvement are described in R], 
we only recall the heavy quark lattice action |q], 



(3) 



(Vq the backward derivative's time component, 
i.e. static quarks propagate only forward in time), 
and the 0(a) improvement term of ^q'*^* in 

{Ar% = At^' + acTii,^n^\(^, + ^*)^h, (4) 

where the 1-loop c^^* = -l/(47r) ^ is used. 

Although in ref. 0] a suitable renormalization 
condition for Aq**** in the SF scheme has already 
been formulated, we slightly modify it in the 
present non-perturbative context for reasons of 
numerical precision. We impose instead 



X{u,alL) = 



/i'-*(L/2) 



_/i x/hh(i/2) 



,1/4 



(5) 



at vanishing quark mass. /I*^^* is a correlator be- 
tween a static-light pseudoscalar boundary source 
and Ag**** in the bulk, /i between two light-quark 
boundary sources, and /{'^ denotes a correlator in 



xa-direction of Wilson lines between static-light 
boundary sources, which is efficiently evaluated 
in MC by multi-hit ||]. Note that the boundary 
wave function renormalizations and a linearly di- 
vergent mass counterterm cancel in eq. (^). 

Fig. |2| displays the continuum SSF after a ^ 
extrapolation of SJv*^' together with an interpolat- 
ing fit. crj^'^* being known, it connects <&(/i) in the 
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Figure 2. Continuum SSF and polynomial fit. 

low-energy regime at fj.ni = l/^m = (2Lmax)~^ 
step-by-step with the perturbative domain, and 
once fj, is large enough for perturbative running to 
set in, we find by numerical integration in eq. (|^): 

*(Ai)/$RGi = 1.088(8) at m = A^m • (6) 

In fig. 1^ we compare the numerically computed 
running of the static axial current in the SF 
scheme with perturbation theory, where also 

enters the analysis. 
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Figure 3. Running matrix element of (^r'^')o in 
the SF sclieme compared to perturbation theory. 

3.1. The total renormalization factor 

Wc still need to relate (A^''*)o(/i) to the 
bare lattice operator. This amounts to com- 
pute Zj^^^ at the (low-energy) matching scale 
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Lm = 1.436 To- Fig. U shows numerical results for 
Z^'^^{go, L/a)\L=L^ and polynomial fits, where 
we also studied the cases cj'^* = and Csw = for 
later purposes. The total renormalization factor 
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Figure 4. Results for Zf'^^ at scale L = 1.436 tq. 
to translate any bare matrix element $barc(5o) of 

/tstat 

^*(.9o) 



Q into the RGI one, $rgi, then reads 

«'RGI 



xZT\g^,L/a) 



(7) 



with <l'RGi/*l'(Aim) a universal part independent of 
the regularization and Z^^^^{gQ,L/a)\L=L^ that 
depends on it. Z$ will be given explicitly in 1^. 

3.2. Wilson data on F-q revisited 

Since we now also have the SF renormaliza- 
tion factor Z|'^**((7o, .^/a)|L=L„, for unimproved 
Wilson fermions (csw = 0) available, it may be 
combined with the universal part <i>RGi/$(/.t) at 
/i = /im and the pcrturbativc conversion factor 

$match(Ai)/'5RGI at SCalc H = ™b,MS ^'^'^^'^ 

to confront it with tadpole-improved estimates 
of the FNAL group |jl^ on the Z-factor, which 
equals <I>matcii(TOb.Ms)/*J'baro(go) in our notation. 
As it turns out (cf . IT] ) , they deviate significantly 
in the relevant range of goj which reveals the im- 
portance of non-perturbative renormalization. 

Just so, we non-perturbatively renormalize ex- 
isting Wilson data 0,0 on $barc = F^'^'^'^J^, 
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w 3/2 ^ _ 3/2 
X rJ $barc = rJ FsF 



and extrapolate them to a = in fig. |^ discarding 
the rightmost point {(3 = 5.7). With ro = 0.5 fm, 
this results in F§^^ = 261(46) MeV at ^ = 'rn^;MS 
in the 'matching' scheme, containing all errors 
apart from quenching. 
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Figure 5. Continuum extrapolation of the non- 
pcrturbatively rcnormalizcd F^^'^° from |I^ , |TI| /. 

4. CONCLUSIONS 

We have performed the scale dependent renor- 
malization of Aq'''* by constructing a non- 
perturbative RG in the SF scheme, and agree- 
ment with perturbation theory at large scales 
is demonstrated. The renormalization factors 
needed to get the associated RG invariants are 
computed with good numerical accuracy 
which is a crucial prerequisite for a controlled de- 
termination of .Fb in the static limit. 

Our continuum extrapolation that uses unim- 
proved data from the literature and also quite 
large lattice spacings leaves much room to im- 
prove on the present result F^^^ = 261(46) MeV. 
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